NONADIABATIC EFFECTS IN A GENERALIZED JAHN-TELLER LATTICE 
MODEL: HEAVY AND LIGHT POLARONS, PAIRING AND METAL-INSULATOR 

TRANSITION 
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The selfconsistent ground state polaron potential of one-dimensional lattice of two-level molecules 
with spinless electrons and two dispersionless phonon modes with linear coupling and quantum 
phonon-assisted (nonadiabatic) transitions between the levels is found anharmonic in phonon dis- 
placements. As a function of these, the potential shows a crossover from two nonequivalent broad 
minima to a single narrow minimum which correspond to the positions of the levels in the ground 
state. Generalized variational approach respecting the mixing of levels (reflection) via a variational 
parameter implies prominent nonadiabatic effects: (i) In the limit of the symmetric E®e Jahn- Teller 
situation they cause transition between the regime of the predominantly one- level "heavy" polaron 
and a "light" polaron oscillating between the levels due to phonon assistance with almost vanishing 
polaron displacement. Vanishing polaron selflocalization implies enhancement of the electron trans- 
fer due to decrease of the "heavy" polaron mass (undressing) at the point of the transition. There 
can occur pairing of "light" polarons due to exchange of virtual phonons. Continuous transition to 
new energy ground state close to the transition from "heavy" polaron phase to "light" (bi)polaron 
phase occurs. In the "heavy" phase, we have found anomalous (anharmonic) enhancements of quan- 
tum fluctuations of the phonon coordinate, conjugated momentum and their product in the ground 
state as functions of the effective coupling which reach their maxima at E(g)e JT symmetry. They 
decrease rapidly to their harmonic values as soon as the "light" phase is stabilized, (ii) Nonadiabatic 
dependence of the polaron mass (Debye- Waller screening) on the optical phonon frequency appears, 
(iii) The contribution of Rabi oscillations to the transfer enhances significantly quantum shift of the 
insulator-metal transition line to higher values of the critical effective electron-phonon coupling sup- 
porting so the metallic phase. In the Eg)e JT case, insulator-metal transition can coincide with the 
transition between the "heavy" and the "light" (bi)polaron phase only at certain (strong) effective 
electron-phonon interaction. 

PACS number(s): 71.38.-k, 42.50.L, 31.30.Gs, 74.20. Mn 



I. INTRODUCTION 

Recently an interest in electron-phonon models 
was renewed owing to high-Tc-superconducting layered 
cuprates and fuUerene compounds which exhibit struc- 
ture instability due to strong Jahn-|Tejl]er effect accom- 
panied by an evident isotope efFectlll^B. Here, nonadi- 
abatic (quantum) fluctuations were found to be impor- 
tant since the majority of the superconducting structures 
exhibited low values of the Fermi energy, comparable to 
the phonon energy and large isotope effect. Nonadiabatic 
fluctuations also appeared relevant in manganese-based 
perovskites at Jahn- Teller distortion assisted by forma- 
tion of JT polaron and causing oxygen isotope effect and 
a coUosal magnetoresistanceO. The mechanism of nonadi- 
abatic pairing was proposed by several, authors : Manini 
et|-|alj, Zheng et alt2l^E3, Kresin et aoB, Pietronero et 
alD for Ceo compounds and recently also for the chal- 
lenging new superconductor MgB^^. Theories of pair- 
ing mechanism based on polarons (bipolarons), including 
also clfiGtipii-electron interactions evoked great interest 
as wellll3"lla. 



Nonadiabatic fluctuations affect the charge transpcpt 
because they reduce the polaron band narrowingt3^t2l 
(i.e. polaron renormalization of the elcctroHjn|aas). Fluc- 
tuations increase near the phase transitionaij^lij destroy- 
ing the phase coherence for weak interactions and shift- 
ing the critical couplings to higher values. However, nu- 
merical sinmlations of diffccejit lattice models .-(Fradkin 
and Hirschlij, Borghi et al.EIl, McKenzie et alta) prove 
a general statement that the resulting quantum fluctua- 
tions are much more pronounced than those |©btained by 
variational -approaches, e.g. bj2-.Zheng et al.t3, Feinberg 
et al.E3, LcEJ and Chen et al.El. Namely, the quantum 
shift of critical values of the electron-phonon coupling 
from numerical simulations to higher values (preferring 
so the metallic phase) was shown to exceed consider- 
ably their variational values. Moreover, the numerical 
simulations of two-level lattice electmn-phonon systems 
in [OJic dimension by Feinberg et al.Ej and by Borghi et 
al.EIl evidence anomalous increase of quantum fluctua- 
tions of the phonon coordinate and of the conjugated 
momentum as well as anomalous increase of their prod- 
uct referring to enhancement of anharmonic effects. The 
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simulations manifest a dramatic increase of phonon fluc- 
tuations and of their product far beyond the standard 
uncertainty principle. The inadequacy of the variational 
approaches was ascribed to insufflcicncy of the squeezed 
coherent phonons (harmonic oscillators) to comprise an- 
harmonic behavior of the phonons in two-level (-band) 
models. l_. 

Since the early beginningpj Jahn- Teller model was in- 
vestigated mainly in chemistry as a prototype model of 
electron-phonon interaction for localized centers in solids. 
It was used to study the instability of the orbitally de- 
generate electronic states of highly symmetric ionic con- 
figuration against ionic |distortions in localized molecu- 
lar centers in crystalscd'Eil. The effect was explained by 
Jahn- Teller theorem about lifting the degeneracy of an 
orbitally degenerate electronic state by symmetry lower- 
ing distortions of nuclear configurations. This effect is 
also involved e.g. in physics of structural phase transi- 
tions in solids doped with Jahn- Teller active ionsE3, for 
mechanisms based on tunneling-between two. electronic 
levels coupled to phonon modescj, in opticala and para- 
magnetic ion spectraBZl. 

Though JT effect was considered as the most represen- 
tative nonadiabatic system, Born-Oppenheimer approxi- 
mation has been extensively used although it is valid only 
in the limit of large local distortions. Inconsistency of the 
adiabatic. acproach for small distortions was analyzed by 
WagneiO'Ea. 

Importance of JT effect in physics was increased due 
to above mentioned discovery of JT effect based struc- 
tural transitions in some high-Tc superconductors and in 
manganese-based perovskites. Therefore, the considera- 
tion is now focused on lattice versions of the JT model. 

In this paper we will investigate an extended lattice 
of JT molecules: site molecules with doubly degenerate 
electron level coupled to two internal phonon modes with 
different symmetries and different coupling constants. 
The antisymmetric mode splits the levels, the symmet- 
ric one couples the levels via phonon assisted on-site and 
intersite electron transitions (Rabi oscillations). As we 
show, these transitions are especially important if the 
difference of the effective potential minima as a function 
of the phonon displacement is of the order of the phonon 
energy (symmetric Jahn- Teller molecules) . We shall use a 
generalized variational approach which will account more 
aspects of nonadiabaticity: Instead of simple squeezed co- 
herent phonons, as a more appropriate variational ansatz 
for phonons we take two-center squeezed coherent wave 
function which accounts for possibility of the phonon me- 
diated coupling of levels. The variational wave function 
of the two-level electron-phonon systems with reflection 
symmetry (antisymmetry) of the coherent phonon states 
related to both the levels (j>+ and </)_ was first introduced 
by Shore and Sandertj as linear combinations + rjcj)^ 
and (jj- +r](j)+, rj being new (reflection) variational param- 
eter. The approach has been widelv exploited and fur- 
ther developed by Wagner et alEJ'EJ for exciton ((211)), 
dimer ((221) and (222)) systems (with the tunneling be- 



tween the levels in contrast to the phonon-assisted elec- 
tron transitions of our model); here the notation (xyz) 
is used for x as number of the levels, y as the number of 
the sites, and z as number of the electrons in a cell). 
The detailed comparison of the ground state energies 
(GSE) of different variational phonon wave fupctions for 
an exciton-phonon or dimer-phonon modelsEJ conflrmed 
the two-peak variational choice as the most suitable one, 
i.e. giving the best fitting of the exact solution in the 
medium and strong coupling regime. 

Similar structure of the wave function was proposed 
for the phonon wave function coupled to two electronic 
states of a double- well potential by Kresin et aM in their 
model of the nonadiabatic origin of the isotope effect in 
high-Tc superconductors. 

In the Section II. we formulate variational approach 
for our model. We find effective ground state potential 
as a function of all (four) phonon variational parameters 
(displacements of two phonon modes, squeezing and re- 
flection (mixing)). 

In the Section III. we compare variational results for 
the ground state energy of the two-center squeezed co- 
herent phonon wave function with the adiabatic ground 
state energy and identify the region of importance of the 
reflection (77) and nonsymmetry (% = P/a) parameters 
to be the region close to x = 1 (E(g)e JT symmetry). 
Special importance of the symmetric (Jahn- Teller) lat- 
tice molecules occurs as a specific condition for which 
the reflection parameter is the most effective. At this 
symmetry, there takes place a nonadiabatic transition 
between the polaron dominantly selflocalized within one 
level, "heavy" polaron of a broad minimum and the al- 
most delocalized, "light" polaron of a narrow minimum 
(with vanishing displacement) oscillating between close 
levels via both onsite and intersite phonon assisted tun- 
neling: Namely, the coherent phonons 2 are accompanied 
by Rabi oscillations of the electron between the levels due 
to the phonon mode 1. The latter virtual phonons medi- 
ate the coupling of polarons which may occupy the levels. 

We compare ground state characteristics of the 
"heavy" and "light" region as functions of pairs of 
competing classical (effective coupling) and quantum 
(phonon frequency, tunnelling) parameters. The self- 
trapping due to electron-phonon coupling competes with 
the lattice transfer of electrons (band width) which it- 
self is being renormalized by the Debye- Waller factor. 
The characteristics of the related insulator-metal tran- 
sition, the value of the gap and consequently the crit- 
ical coupling are determined by the complex interplay 
of the transfer supporting the metallic phase and of 
the electron-phonon coupling supporting the insulating 
phase. The shift of the critical line to higher values of 
the critical couplings due to quantum effects is discussed 
in the Section IV. 

In the Section V. we investigate anomalous behavior 
of the squeezed ground state quantum fluctuations of 
the canonically conjugated phonon coordinates and their 
product on the effective coupling (/i) and effective po- 
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tential asymmetry (x) parameters, namely the strong an- 
harmonic fluctuations which reach their maximum values 
again a,t x — 1 • They decrease to the harmonic oscillator 
values for X > 1. A related model is a lattice of two-level 
dimers with one spin electron at each site ((222) lattice 
model! studied by exact numerical methods by Borghi 
et alO. The present model differs from that one, but 
concerning the quantum fluctuations of the phonon co- 
ordinate and momenta qualitatively the same results as 
for their model are expected. 

We remark that variational methods are widely used 
for electron-phonon lattice models. 

Fully analytic nonvariational approach with noncon- 
servation of the number of phonons was pjerformed for 
a local {D = 0) dimer by Weber-MilbrodtE3. Unfortu- 
nately, in the extended (D=l) lattice model this method 
acquires extreme mathematical-complexness. It was used 
as a basis for numerical studyO of the above mentioned 
lattice model. 



II. EXTENDED (LATTICE) GENERALIZED 
JAHN-TELLER MODEL 

We investigate ID lattice of spinless double degener- 
ated electron states linearly coupled to two intramolecu- 
lar phonon modes described by Hamiltonian 

H = ^ iblhn + l/2) + Y.(a{bl+h,n)'yzn 

n,i=l,2 n 

-PibL+b2^n)a,n)-J {Rl,l + R-l.j)In. (1) 



a 



n,j=l,2 



where bi^nji 



1, 2 are phonon annihilation oper- 
ators, and the Pauli matrices (T;„ represent two- level 
electron system. They satisfy identities [o';„,crj„] — 
io'kn, I — x,y,z, representing 1/2-pseudo-spins related 
to the electron densities in a usual way, i.e. Uxn — 

2(''l:n''2,n + C2 „Ci^„), fJyn = (''l,n'^2,n — ^j ^Ci^n), 
^zn = 2('^l,n''l." ~ '^2,n''2,n)i In = 2('^l,n''l." '^2,n''2,n) 

is a unit matrix, and Cj^n are electron annihilation oper- 
ators. The operator R±ij = e^'^" of the displacement 
by a lattice constant ±a acts in both the electron and 
phonon space, R±i,jfn = /n±i-R±ij- 

In terms of the creation-annihilation electron and 
phonon operators the Hamiltonian can be cast as follows: 

^^Y^^Y ^^lnb^n + J) + ^{um ^ «2„)(6L + 6i„) 



- 2 ('^1«'=2« + 4«ClrO(4« + b2n) 
3 = 1,2 

For f3 — ~a, the interaction part of (^ 



bin + bi 

'An + b2n, 



bin + b2n 

-(bL + bm) 



(3) 



yields the rotationally symmetric E ® e for nJi with a 
pair (an antisymmetric and a symmetric under reflec- 
tion) of double degenerated vibrations. This is, e.g., 
ions with cP configurations in high- Tc 



the casCfOfrJ 
cupratesEJ'O. 

Taking a ^ (3 removes the degeneration of the vibronic 
states breaking the rotational symmetry of the electron- 
phonon interactiona—tlie model still staying within the 
class of JT modelsBEa-B. 

The dispersionless optical phonon mode bi splits the 
degenerated unperturbed electron level {j = 1,2) while 
the mode 62 mediates the electron transitions between 
the levels. This latter term represents phonon-assisted 
tunneling, a mechanism of the nonclassical (nonadia- 
batic) nature as well as is the pure tunneling in related 
exciton and dimer models. 

Evidently, Hamiltonian (Q) {a 7^ [3) is reflection- 
symmetric, G'^^'^G'^P'^^H = H, 

|2) = G(^')|1), (G('='))2 = 1, 

Gt'\bl±b,,,) = -{bl±bu.)Gt\ {Gt''r = l. (4) 

where G^^'*'* = exp(iiTb\^bin) is the phonon reflection op- 
erator. While the phonon 1 is antisymmetric under the 
reflection, phonon 2 remains symmetric. 

In addition, the transfer part of (^) exhibits SU{2) 
symmetry of the left- and right-moving electrons (holes). 

Let us note that the quantum phonon assistance of the 
electron tunneling (/?-term in (||) and (|l|)) constitutes the 
difference of the model from the related dimer and exci- 
ton quantum models where instead of /? X](4n + b2n)o'xn 

n 

where A is the distance be- 

291,1 ■ ^" 



of (|^) there stajAds-.Aax 
tween the levelsE^'tJ. 

The local part of (|l|) can be diagonalized in electran. 
subspace by the Fulton-Gouterman unitary operatorL^ 

Un = t/2,nC^l,ri, whcrC 

U^.n = ^ (i\'_g;"J ' G.,n = exp(z^4&„) EE G?„'\ 

(5) 

as follows 

HL=J2^nHonU-'^n Y (blhn + l) 
n n,i=l,2 ^ ^ 

+ 1 Y^^ibln + bi,n) - /3(6L + 62„)Gi^„)]/„. (6) 



On the other hand, in the transfer term 



(2) 



Ht^-jJ2 iVn,lRl + Vn,-lR-l) 



(7) 



there appears a nondiagonality 
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Vn,±l — [(1 + Gi„Gi„±i)/„ + (1 — GinGi^n±i)G2nO'zn] 

X [(1 + G2nG2n±l)In + (1 ^ G2nG2,n±l)'^xn] ■ (8) 



Here, Pauli matrices transform as UicrxU: 



GiCTz 



U^a.U-^ = a,, and C/,(6| + h)U-^ = {bl + h)2a^. 

The diagonal terms of (^) represent the polaron trans- 
fer within one level while the off-diagonal ones represent 
the interlevel polaron transfer through the lattice. Evi- 
dently, the contribution of the off-diagonal terms propor- 
tional to 1 — GinGi^n+i is much smaller when compared 
with those proportional to 1 -I- GmGi^n+i- 

Because of nonconservation of the number of coherent 
phonons, they are able even in the ground state to assist 
electron transitions between the levels. In the Hamilto- 
nian ^, the operator Gi„ ~ (— (||), highly non- 
linear in the phonon-1 appears mediated by phonons 2. 
It introduces multiple electron oscillations between the 
split levels mediated by continuous virtual absorption and 
emission of the phonons 1. The effect is analogous j-to 
Rabi oscillations in quantum optics due to photonstj. 
Let us note that Rabi oscillations assist both the inter- 
level onsite and intersite electron transitions mediated by 
the electron transfer T. 

For various local electron-phonon models, a number 
of phonon variational wave functions using coherent and 
squeezed coherent phonona A¥as_nmposed in order to 
reach the best ground stateB^EilMEl. r-| 

Specially, for two-level systems Shore et al.E^I proposed 
variational eigenfunctions in the form of symmetric and 
antisymmetric combinations of the reflection symmetric 
components related to both the levels. Each of these com- 
ponents was chosen in a form which accounted for mixing 
with the reflected state via new variational parameter. 
This choice was stated to be the best variational wave 
function, i.e. yielding the lowest ground state energjt. 
The method was further developed by Wagner et al.E3 
for dimers and excitons. We shall adopt this approach in 
what follows by taking 



1 

1 



C 



. (-01 ± V'>P2) ± (/>- (-02 ± VPl)] , 



(9) 



with the upper sign for the ground state $i and lower 
sign for the excited state $2. 

In (^), = <j)+ +ri4>- and 02 ~ 4>- +'70+ are phonon 
wave functions related to two levels with mixing of the 
bare 0+ (0-) and reflected 0_ (0+) parts mediated by 
the variational parameter 77. 

Further, ij^i and "02 are components of the electron 
state vector related to the upper and lower level, respec- 
tively, 



= cl„|Oei), V2N = 4n|0ez) 



(10) 



where jOez) is the electron vacuum. Note that the last line 
of (H) allows us to interpret alternatively the parameter 
rj as reflection parameter of the electron states as well. 
The squeezed coherent phonon wave functions 



0± [llrn{n),^2ni{n),rrn{n)] 

Di,±[lim{n)]D2 [72m("-)]'5'i [rmin)]\Qph), 



(11) 



are related to lower (0+) and higher (0-) level; \Qph) is 
the phonon vacuum. 

The generators of the coherent phonons arc functionals 
of related displacements 

^i,±[7im(")] = exp^±[7i„(n)6|,„ - 7i*m(")&im], (12) 

m 

D2[l2m{n)] = exp^[72™(n)4„j - 72m(»^)^2m], (13) 
m 

as well as the generator of squeezing 

5i[r„(n)] = exp^r„(n)(C - b\J. (14) 

m 

Respective generator of the mode 2, S2[r2m{n)] — 1, be- 
cause Hamiltonian (H) is linear in 62, b\ and r2m{n) = 0. 
The normalization condition then reads 



where 



and 



Ci„(n) = l-hry^-h2?7exp(-2|7i„(n)p). 



limin) = 7im(n)exp(-2r„i(n)). 



(15) 



(16) 



(17) 



Generally, ^jm{n) is a nonlocal quantity which rep- 
resents the displacement of the mode j at the site 
m due to an electron at the site n, ljm{n) = 
X] Ijqi^) exp{—iqma). We take 



e:xjp{iqna) -> 7jm("-) = IjSr, 



(18) 



where 7^ are independent of n. Eq. (|1^) indicates that 
the phonon displacement accompanies the electron at the 
site n. The same is valid for the squeezing, rm{n) = 

Full Bloch solution to the transformed lattice Hamil- 
tonian H = Hl + Ht, (H), (0) is chosgHpas a generalized 
Fulton-Gouterman variational ansatzBE3 in the vector 
form (Eq. (|)) 



J.i:„p(.:.„e,(^(»)*W), m 



where the electron parts ipi are defined by ( p^ ) and the 
phonon parts 0i(n) by (^ and (ll)-(p^. 
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We have investigated the model (||) in a former paperES 
with respect to the stabiUty of a sohton ground state 
against quantum fluctuations. However, the variational 
wave function we used did not account for the two-center 
nature of the wave function here respected by the reflec- 
tion parameter 77 (^. 

In what follows the ground state energy as function 
of optimized variational parameters rj, 71,72,^ will be 
determined. 



III. THE GROUND STATE 

When considering the ground state oi H = Hl + Ht, 
and (0), we omit the electron and phonon dynamic 
terms by setting fc = 0, 5 = 0. 

By averaging H = Hl + Ht ( fl), ( [t])) over the phonon 
wave functions with (^) and (|ll|) one obtains for the 
site Hamiltonianl^) (Appendix B) 



(20) 



where 



Hi, 



Hph^^{cos\i^r + l)I (21) 

l + ^2_2^^-8..exp(-2|7iP) 2, , 0,., ,2, ...^ 
l + ,2 + 2,ycxp(-2|7iP) ' + ^''^ 



Hn + Hp, 



l + 772 + 27;exp(-2|7i|2) 



/3[(l + 772)exp(^2|7i,|2) + 27?] 

l + r72 + 2ryexp(-2|7iP) , I 

where iJa and iJ/j are a- and /3-dependent parts of the 
interaction i?i„f (|2^). From the transfer Hamiltonian //j- 
(^ there remains 



-rexp(-Ty)M 



2{HtiI + HT2i^y + HxiCFz + Hncfx) 



(24) 



The Debye-Waller factor exp(— W)Af in ( |2^ ) is given 
in the Appendix B. The differential part of the Debye- 
Waller factor W in the ground state is zero. For M it 
holds 



(25) 



M = [(1 + El)I + (1 - El)E2<J,,,] 
X [{l + El)I+{l-El)<j,] , 

where Ei are given by (Appendix B) 

El = exp (-2|7i|2) , E2 = exp (-2|72p) . (26) 

Expressions Ei ( p6| ) in the ground state are independent 
of q and n because of the form of jiq{n) ~ 7^ exp(— igna) 



(|l8|). This substantially simplifies subsequent calcula- 
tions leaving ii^i and Ci as functions of independent 
on n. The transfer terms Hti, i = 1,2,3,4 in (|2j) with 
( ^ ) are expressed as 



H 



T 



Tl 



{1 + El){l + Ei), 



T 



Ht2 ^ - Et){l - Ei)E2, 

HT3 = -^il-El)il + El)E2, 
HT4 = -^{l + Ef){l-El). 



(27) 



Here, Hti,Ht3 and Ht2,Hta are diagonal and off- 
diagonal terms of (U), respectively. 

The effective polaron potential in ^ (with (|2^)-(|26|)) 
as a highly nonlinear function of 71 and 72 is visualized 
in Fig. 1. 

The potential exhibits two sets of minima related to 
two competing ground states: 

(i) Two nonequivalent broad minima related to both the 
levels (^2]) at ±71 ^ and 72 close to 0; 

(ii) One narrow minimum at 71 close to 0, 72 > 0, where 
both levels approach close together. This minimum de- 
velops at growing /3; evidently its behaviour depends also 
on T value because of nonlinearity of the Debye-Waller 
factor M and (|6|)). 

The ground state energies related to these two sets 
of competing minima were calculated numerically. The 
result of the numerical minimalization of the diagonalized 
form of energy (EOl) 



E — Hph + Hp -f Ht^ — [{Ha + Hts 
can be written symbolically as 



H 



2 1I/2 



(28) 



EG{^J■, X, ^/T) = E{{'^ig, J2q,r, 77}mm)/T- 



(29) 



Here, the index min denotes the optimized values. The 
model parameters 



(30) 



are parameters of the effective interaction, asymmetry 
and nonadiabaticity, respectively. Energy Eq (^9|) is in 
scaled units, renormalized by T. The results of the nu- 
merical evaluations of the ground state energy ( p9| ) are 
depicted in Figs. 2a, 2b different model parameters. 

One can distinguish there two regions depending on 
the value of x with different behaviour of ground state in 
each of them: 

(i) The ground state pertaining to the lower broad mini- 
mum at 7i < and 72 ~ with a small reflection part at 
71 > is referred to as a "heavy" region. It corresponds 
to predominantly intralevel " heavy" polaron which is rep- 
resented by a two-peak wave function, both peaks rep- 
resenting a harmonic oscillator exp[— (x =F 7i)^]) dis- 
placed by the value ±71 (as it is seen from the form of 
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the ansatz (|12|)). The "heavy" polaron is dominant at 
X < 1 where the broad minimum 71 < is dominating, 
(ii) Close to X = the energies of two minima go very 
close together (their difference is of the order of the 
phonon energy Q). They drop to one narrow minimum 
which represents a new ground state. Close to x — ^, 
continuous transition to a new ground state occurs. It 
stabilizes at x > li where (optimized) 71 value is close 
to and 72 > 0. This region is referred to as a "light" 
polaron region because, owing to the abrupt decrease of 
7i « at X > Ij the effective mass of the intralevel 
polaron drops to almost its free-electron value, i.e., the 
polaron selflocalization vanishes in the "light" region. Be- 
cause of this " undressing" , the transport characteristics 
of the excited electron would increase. Moreover, due to 
the tiny distance between the levels, their coupling takes 
place by the exchange of virtual phonons 1. At suitable 
conditions, in the excited state, when both levels are oc- 
cupied by electrons of opposite spins, the mechanism of 
virtual phonon exchange implies the pairing of electrons, 
i.e. formation of "light" bipolarons . 

The ground state energy, especially its behavior depen- 
dent on pairs of parameters X: ^ smd x, fJ- is illustrated 
in details in Fig. 2. While being weakly fi-dependent, 
the energy strongly decreases with x inside the "light" 
phase (Fig. 2a). The position of the phase line is slightly 
shifted from x = 1 at 51 = to higher values of x with 
increasing f2. This is consistent with the fact that the 
phonon fluctuations are most effective when the differ- 
ence of the energies of the phases is of the order of the 
phonon energy. 

The ground state energy in the "heavy" region (x < 
1) is independent of Xj its decrease inside the "light" 
region (x > 1) is dependent on the effective coupling fi 
(Fig. 2b). The energy decrease due to fi is stronger in 
the "light" phase (depending on x) than in the "heavy" 
phase. 

The region of importance of reflection measure rj and 
r is illustrated in the Fig. 3. There we show the dif- 
ference between the ground state energy with r and rj 
omitted and ground state with four variational param- 
eters calculated numerically. However, contribution due 
to r Ecirj, r = 0) — Eg ~ 2.10"'^ of the maximum value in 
Fig. 3. The said parameters turn out to be most relevant 
in E(g)e JT case (x = !)■ 

The effects of rj and r on the ground state is apparent 
only at moderate couplings in the "heavy" region, reach- 
ing their maximum at x = 1- The narrow minimum (the 
"light" region) is resistant against 77. 

In order to demonstrate important properties of dis- 
placements (/i and Q dependence) it is sufficient, except 
of the region of importance of the reflection (Fig. 3) 
inside the "heavy" region apparent for 72 (Fig. 4), to 
calculate them in the limit 77 = 0. We shall use Eq. ( |2^ ) 
approximated for the cases of both " heavy" and " light" 
polaron and obtain implicit expressions for 7^ which how- 
ever are good to visualize their behaviour in both regions: 

(i) "heavy" polaron (72 small, E2 ~ 1,Ht2,Ht4 ~ 0): 



71 



-ae 



4r 



2TEf 



(31) 



72 is small except of the region of fluctuations visualized 
in Fig. 4. One can see, that the selflocalization due to 
phonons 2 in the "heavy" region — 7I clearly implies 
the " cave" in the ground state energy due to the reflection 
measure 77 (anharmonicity) of the ground state (Fig. 3). 
(ii) "light" polaron (71 smafl, Ei ss I^Ht^^^Ht^ « 0): 



-ae 



4r 



71 



72 



2TEl 



2nT 



10" 



f3 



2TEl{2- 



(32) 



In the "light" region the fluctuations of 71 are missing 
as well as the fluctuations of the energy. This is consis- 
tent with the above result of the resistance of the narrow 
minimum against 77. 

For both " heavy" and " light" polaron a dependence of 
7i on the nonadiabaticity parameter f2/T appears. It im- 
plies the dependence of the Debye- Waller factor and con- 
sequently of the polaron mass on the phonon frequency 
fi. This can be thought of as an analogy of the isotope 
effect at zero temperature. 

In the "heavy" phase, electron transitions mediated by 
phonons 2 to the upper level enhance fluctuations ~ 72 
which mix phonons 2 with phonons 1 and contribute to 
the fluctuations of the ground state energy of the heavy 
region. This is the reason for the similarity of the results 
in Figs. 3 and 4. 



IV. METAL-INSULATOR TRANSITION 

To illustrate the effect of various interactions let us 
consider first a simple two level case (T = 0). The ground 
state will be obtained by the numerical minimalization of 
energy 



E^ — Hph + Ha + 



01 



(33) 



where Ha and Hp are a and (3 dependent components of 
the interaction term Hint (M3|). The result is of the form 



Eg = ^(cosh4r0 + 1) - Hpa - Hpp, 



(34) 



where the label specifies the optimized variational val- 
ues. In (p^), Hpa and Hpp are contributions of the selflo- 
calization energies of two coupled polarons proportional 
to /2Q, and 0^/2^, respectively. The splitting of the 
levels is implied to be 



A" 



(35) 



Respectively a lattice of JT molecules (T ^ 0) shows 
either metallic or insulating phase behavior depending 
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on whether the spUtted bands overlap or not. Usually 
the resulting situation depends strongly on the phonon 
renormalization of the gap and of the electron band in- 
cluding phonon fluctuations. 

Our aim here is to investigate the metal-insulator tran- 
sition due to the competition of the polaron localization 
and the delocalization due to the transfer including the 
phonon assisted electron transfer effects. Namely, the 
transport terms Hti compete the level splitting (35) be- 
cause of the complex interplay of the diagonal and off- 
diagonal contributions to the band width. In order to find 
the value of the gap as an order parameter of the metal- 
insulator phase transition it is necessary to perform usual 
renormalization: SU{2) symmetry of the right- and left- 
moving electrons and holes of the transfer part of Hamil- 
tonian (^ makes it possible to split electron operators 
onto an electron and a hole part, 



Ci,-, i = 1,2, 



(36) 



the signs ± indicating directions of the electron motion. 
In view of the reflection symmetry of the two-level local 
part of Hamiltonian (|^) against the level E = in the 
middle between the levels electron operators (^6|) can be 
rewritten in a renormalized form 



Cl 



~ ^1,+ ''1,+ ' 



C2 



(37) 



This renormalization implies the change of the sign of T 
(Ht,) for holes (£; < 0). 

There are two pairs of electrons and holes, cf''_^, 



and Cj*^^, c'l'^l! , i.e., two state vectors ( 



and 



-<J- 

n 



which satisfy equations 



TT I f+ 

.9- 



-A/2, 0\ f-g*_ 
0, A/2 J [ U 



Hn 



n 



A/2, 
0, -A/2 



9- 



= 0, 



(38) 



(39) 



where Hm results from the minimalization of ( pO| ) with 
respect to the variational parameters, (By), 



Hm — 



-Eg, 
0, Ea 



(40) 



and Eg given by (|2£ 

The solution for the gap A reads finally (see Figs. 2b- 
4b) 



A = 2|£;G(T^-r)| 



(41) 



independent of x up to x = 1- In the "light" phase, the 
critical coupling ^L{fl,x) < ^J-ni^) decreases with in- 
creasing X- The point (x = 1, ^J■Hi^)) is the only point 
of coincidence of the "heavy" -"light" transition and the 
I-M transition. At small /i the fluctuations support the 
metallic phase in the "light" region (Fig. 5a), while at 
small X (and sufficiently large /i that the insulating phase 
could be stable) they support the metallic phase in the 
heavy region (Figs. 5a, b). 

The I-M transition can appear in the " heavy" region at 
large fi ( |30| ) and small (Fig. 2a). From Fig. 3 (medium 
values oijT) one can locate the region of relevant reflection 
r] into the metallic phase, where the bands overlap. The 
"light" region (x > 1) is resistant against ry: no similar 
effect as that one shown in Fig. 3 there exists. There- 
fore, in the "heavy" region the metallic phase is strongly 
supported due to the fluctuations at weak couplings /i. 
"Heavy" phase is therefore metallic in the broad range 
of parameters except of strong couplings (Fig. 2b). At 
certain (large) fj, the I-M transition drops close to the fine 
X = 1, so that the transition coincides with the "heavy" 
to "light" polaron transition. With decreasing /i the I- 
M transition line moves to larger x (Fig- 5b): the M-I 
transition moves away from the "heavy" -"light" transi- 
tion into the "light" region. In the "light" region the I-M 
transition line moves to smaller /i when compared with 
the heavy region (Fig. 6). 

As one can also see from Figs. 6 and 8, the I-M 
transition line in the phase diagram /i, f2 is strongly fl- 
dependent. The phase line in the plane (/i, f2) between 
the phases shows the broadening of the metallic phase 
with growing Q, i.e. the metallic phase is supported by 
the quantum fluctuations. Respective ground state en- 
ergy slightly increases with Q. 

To make effects of different contributions more 
transparent, we can first calculate analytically effect of 
the diagonal terms Ht^ and iJ^s ■ (Let us note that ne- 
glecting the off-diagonal transfer, i.e. Ht^, i = 2,4 in 
(P7|), means neglecting the additional splitting of the 
bands. More precisely, it means neglecting the off- 
diagonal terms in Eq. (|^) ). 

Minimalization of the energy (pq) with = Ht^ = 



E' 



ph 



H' 



(42) 



implies the ground state of a similar structure as that 
one of (ph. 



Eh 



(cosh4r^ + 1) 



(43) 



and implies a condition for the stability of either the insu- 
lating phase, if Eg{T —T) > 0, or the metallic phase 
in the opposite case, where A = 0. 

Fig. 5. depicts both I-M and "heavy" -"light" transi- 
tions. In the "heavy" phase, at small fl the I-M tran- 
sition line is determined by the critical coupling /i//(f2) 



where the label 1 specifies the optimized variational val- 
ues. Respective splitting yields 



Ai = 2\Eh - 2{H}r, 



(44) 
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The qualitative picture of the role of different contri- 
butions in the formation of the gap is given in Fig. 7. 

The exact phase diagram is compared with the result 
of (^J) in Fig. 8. Here one can see that the main con- 
tribution to the shift of the phase diagram in support of 
the metallic phase is on account of the off-diagonal elec- 
tron transfer terms Ht^ and Ht^ related to phonon 
assisted Rabi oscillations. These contributions play im- 
portant role in the "light" phase (Fig. 2b). 



V. SQUEEZED QUANTUM GROUND STATE 
(VACUUM) FLUCTUATIONS OF PHONONS 

In Fig. 9 we display numerical results for the coordi- 
nate and momentum fluctuations of the phonon in the 
"heavy" region (x < 1) of the ground state. Evident 
anharmonicity peaks in the region of relevance of the re- 
flection (Fig. 10) are due to the mixing with the phonon 
of the higher level. The peaks are related to the ground 
state (squeezed vacuum) fluctuations of the phonon co- 
ordinate Qi — ^{b\+bi) 



■ exp(4r) 



47? - 



4(1- 



(l + 7y2 4.2,yexp(-27f))2 



(45) 



7i defined by (p7|). For the variational parameters ji,ri,r 
we have to take their optimized (ground state) values. 
For the related momentum Pi = --^{b\ — bi), one obtains 
the ground state fluctuation 



;Pi)2 = -exp(-4r) 



(46) 



where again the optimized value for r is meant. 

Their product A^QiA^Pi for the oscillator in the 
ground state (Fig. 8) exhibiting a peak in the region of 
relevant reflection rj is the evidence for its anharmonic- 
ity due to the coupling with the higher level oscillator. 
The reflection variational parameter rj is efficient only in 
the "heavy" region (x < 1), for x — ^ the fluctuations 
decrease close towards the classical value 0.25 of the har- 
monic oscillator near the narrow minimum: the stable 
ground state of the "light" polaron is achieved at x = 1- 
In the "light" region the oscillator remains close to har- 
monic one; no analogous effects to those in the "heavy" 
region occur. The results displayed on Fig. 8 stay in ev- 
ident correlation with the region of relevance of the re- 
flection parameter ry of the Fig. 3 (see its detail in Fig. 
10 for the corresponding region of parameters). Similar 
effects of anharmonicity were found by nuoaerical simu- 
lations in the related model by Borghi et al.tZl mentioned 
in the Introduction. 



VI. CONCLUSIONS 

The selfconsistent polaron potential provided by the 
Holstein intralevel and interlevel electron-phonon cou- 
plings in our two-level lattice model is highly nonlinear 
function of phonon displacements (Fig. 1). As a conse- 
quence, there occurs a competition between the regime 
of two nonequivalent broad minima at ±71 and 72 close 
to related to two electron levels and the regime of one 
narrow minimum at 72 > and 71 close to when both 
broad minima collapse to a single one, so that both levels 
drop to that minimum. The broad minima dominate at 
X < 1, the narrow one at x > 1- At E®e Jahn- Teller 
symmetry (x = 1) energies of the broad and narrow min- 
ima coincide and the phonon-assisted tunneling due to 
the nonadiabatic fluctuations reaches its maximum close 
to this limit. When x is approaching to 1, there occurs 
pairing of the levels due to the interlevel onsite and in- 
tersite polaron tunneling mediated by the exchange of 
virtual phonons. The narrow minimum is suppressed if 
the electron transfer energy T decreases. 

Hence, we have identified two regions of stability: (i) 
the region of the dominating "heavy" polaron (x < 1) 
related to the broad minimum at negative 71 (the ab- 
solute minimum) and 72 « and (ii) the region of the 
dominating "light" polaron (x > 1) related to the nar- 
row minimum at 71 ~ (Fig. 1). The almost vanishing 
polaron selflocalization in the "light region" ("undress- 
ing" ) means enhancement of the polaron transfer due to 
respective decrease of the polaron effective mass. Even 
in the "heavy" phase the effect can be present due to 
fluctuations of the "light" phase (Figs. 3 and 4). 

Close to the transition from " heavy" to " light" polaron 
at X = Ij a continuous transition to new ground state that 
is stable at x > ^ occurs (Fig. 2). In this ground state 
the formation of "light" bipolarons is possible which are 
characterized by the enhancement of the effective transfer 
due to the vanishing selflocalization. The pairing of the 
polarons is caused by virtual exchange of phonons 1 me- 
diated by phonons 2 between the electrons in both levels 
at suitable configuration. This effect might be of interest 
from the point ofr-wiew of bipolaron mechanisms of the 
superconductivitytil. The TM transition and " heavy" - 
"light" phase coincide in the crossing point of the lines 
X = 1 and /i^r(il) (critical line of I-M transition in the 
"heavy" phase. Fig. 5a). 

The Debye- Waller factor (^5|) shows nonadiabatic de- 
pendence on the phonon frequency through 7^, (|3l|), 
(^). Then, the effective transfer parameter of electrons 
is respectively increased and the electron effective mass 
decreased. This effect reminds on the "isotope effect", 
however, at zero temperature. 

The metal-insulator transition occurs in both regions: 
"heavy" region is the region of dominating metallic 
phase, where the transition to the insulating phase can 
occur only at sufficiently large /j, (Fig. 5). On the other 
hand, in the "light" region the I-M transition can occur 
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at relatively small fi (Figs. 5 and 6). For certain fi the I- 
M transition can coincide with the "heavy-light" polaron 
transition at x = 1 (Fig. 5a) . The interlevel and the com- 
bined interlevel and intersite electron Rabi oscillations in 
the ground state, are identified to be the main reason for 
the shift of the I-M transition line to stronger effective 
couplings, i.e. in the support of the metallic phase (Figs. 
5 and 7). Fig. 5a illustrates the existence of the "light" 
polaron in both the insulating and metallic phase. 

Fluctuations of the phonon conjugated coordinates and 
their product exhibit peaks (Fig. 9) in the region of pa- 
rameters where the effect of rj (Figs. 3, 4) is relevant. 
Consequently, the effect appears only in the "heavy" re- 
gion reaching its maximum at x = 1 and growing with 
fi. The peaks evidence for the strongly anharmonic be- 
havior of the " heavy" polaron in the region of relevance 
of rj. In the "light" region, there occurs a degeneracy 
over T] (no effect of rj on the minimum of the potential), 
i.e. no analogous effect of the fluctuations, but nearly 
a harmonic oscillator is observed. This is obviously due 
to the specific form of the variational ansatz, which is 
inspired by the shape of the "pure" heavy polaron (i.e. 
with X << !)■ Therefore the ansatz (||), reflecting the 
essentials of the "heavy" polaron, is proper for describing 
the heavy region, but when concerning the " light" region, 
it presents, in some extent, a heuristic extrapolation. We 
could have started from another ansatz underlying the es- 
sentials of "light" polaron (which could be, for example, 
inspired for the "pure" case of x >> 1) in order to in- 
vestigate the "light" region and gain some extrapolation 
towards smaller x- Naturally, also the transition in E(g)e 
Jahn- Teller case (x — 1), needs further considerations. 
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Appendix A. 

We shall use following formulaSS 

D{j)S{r) = S{r)D{j), 7 = je~^^ (Al) 



= exp(^--[7™(n + l)-7(^)]'j , (^3) 
where 7 = 76^^''. 

For averages with using virtual Fock states \n) we use 
formulas 

{0\D-'S-\-lf''SD\0)^,{^\i~lf''\^), 



(A4) 



r!=0 



where |7)<, = SD\0) = 5|7) and. 



(n|7) 



1 / \ "/2 

1 I 7 



' (n!/i)i/2 \^2^ 



X exp 



2 2^l 



(A5) 



Here, /i = cosh2r, v = sinh2r, fi^ — = 1, here v is 
real. 

Appendix B. 

The average of the Hamiltonian in (^ over the states 
( |l9| ) referred to one site results in 

~ Hph + Hint , where 



N 



Hph = 2CN'^ ^ Ci„,„i(cosh(4r™(n)) + l)/„ 



m.n 



(lV-2e-«'^'"("),7exp(-2|7i„(n)n)|7i„(n)p/„ 



n 



Hint = Ha+ Hfj 



Cl (")|72m('T-)| In, 



1 v-^ r 



m 



CN 



2 2^|"(1 ~ )llni<7zn 



_0 

2 

and 



[(1 + 77^) exp(-2|7i„|2) + 277]72™/„ } H.c., {B2) 



S~\r)bS{r) = b cosh 2r + b^ sinh 2/ 



{A2) 



Ht = --^ECi„(n)exp(-iy„(n))M„(f 



(0|5-i(r„(n))D-i(7,„(n))2?„(n + l)^(r,„(n + 1)) 



= 2{HtiI + HT2i(Jy + HT3<y, + Hta<Jx). (B3) 
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Calculation of the Debye-Waller factor exp(— W)M 
(H) and of (H) is based on the formulas (A3)-(A5). We 
approximated the average of a product in the transfer 
term by a product of averages. In the continuum limit, 
from (A3) we receive for W{n) in (B3) 



2N ■ 



dn 



d-y2m{n) 



dn 



m 



With the help of (A4) and ( A5) one obtains for iJ,,, (n) 

,m+l Gi 771/ 7 

Ern{n) = (0„i+i,0„,.|5"^(r,„+i(r7,))5"^(r,„(r7,)) 



xl?-i(7i,„,+i(n))i?-i(7i„(n))(-l)^i-+i''^-+^ 



x(-l)''i'"''^'"^(7i.m(n+l))-D(7i,™+i(ri+l))^(r™(n+l)) 



S{rm+i{n. + l))|0m,0m+i) = exp ^-^ [|7im(n)|^ + 



+ |7lm(" + 1)P + l7l,m+lWP + l7l,m+l(»^ + 1)1^ 



+7lm(^^ + l)7lm('^) + 7l,m+l(?^ + l)7i,m+i(n) 



+7lm(")7l,m+l("-) + (n + l)7im('^ + l)]) 

In the continuum limit (B5) yields Em{n) 
exp(-2|7„i(n)|2). 
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VIII. FIGURE CAPTIONS 

Fig. 1. Potential for x = P/a = 0.5, = 2 and 
n = 0.05. 

Fig. 2. The ground state energy ( p9[ ) in the plane 
ai ^ = 2.5 (a) and (b) in the x-M plane at = 0.5. 
The "light" polaron is evidently much more sensitive to 
effects of X and /i than the "heavy" polaron. 

Fig. 3. The range of relevance of the reflection param- 
eter 77: Difference A£'g = -Eg(?7 = 0,r = 0) — E'c, for 
X = 1, Eg is the exact ground state (|9|). For x < 1, 
the difference increases with x reaching its maximum at 
X = 1. At X > 1 it drops to zero, i.e. the narrow mini- 
mum "light" phase is resistant against r/. 

Fig. 4. Displacement of phonons 2 in the "heavy" 
region 7I; it evidently causes the "cave" due to reflection 
■q of the ground state energy in the "heavy" region (Fig. 
3) , 

Fig. 5. (a) The gap in the plane of /i and x for 
= 0.5. While the critical coupling fini^) for I-M tran- 
sition in the "heavy" phase is very slightly dependent 
of X, ^J'L{^^x) in the "light" phase is x dependent, (b) 
The shift of the metal-insulator transition in the " heavy" 
region due to quantum fluctuations at strong coupling 

Fig. 6. The gap in the "light" region as function of H. 
and fi for x = 1-5. The fluctuations increase with Q. 

Fig. 7. The competition of the level splitting and band 
width T terms (T = TM (||)). Aq is splitting of the 
levels (HJ), Ai neglects the off-dia gon al terms Ht2,Ht4^, 
(^). A is the exact gap given by (Hll)- 

Fig. 8. Comparison of the phase diagram of the ex- 
act model (^ ) (solid line) with that one neglecting off- 
diagonal transfer (dashed line). The metallic phase oc- 
curs on the left and the insulator on the right of the 
respective line. 

Fig. 9. Quantum fluctuations of (a) the coordi- 
nate A^Qi, (b) momentum A^Pi and (c) their product 
A^QiA^Pi. (d) Detail of Fig. 3 underlying the region of 
parameters of Figs, (a)-(c). 
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